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An analytical model of the charging of a spacecraft emitting electron and 
Ion beams has been applied to the case of large spacecraft. In this model. 
Ionization occurs In the sheath due to the return current. Charge neutraliza- 
tion of spherical space charge flow Is examined by solving analytical equations 
numerically. Parametric studies of potential of large spacecraft are performed. 
As In the case of small. spacecraft, the Ions created In the sheath, by the 
returning current play a large role In determining spacecraft potential. 


INTRODUCTION 

The potential difference created between spacecraft-ground, and. the ambient 
plasma during the ejection of. a beam of electrons from a sounding rocket pay- 
load. In the Ionosphere (ref. 1) has. been found to be much less than had origi- 
nally been theoretically predicted (ref. 2). To determine the reasons for 
this limited potential difference, large-vacuum-chamber tests were conducted 
In which electron and Ion currents were ejected from a payload Into a simulated 
Ionosphere. 

As a plausible explanation to the observed current voltage behavior, 
sheath Ionization models. (refs. 3,4) for small spacecraft have been studied. 
Mhen an electron beam Is emitted from a spacecraft, ambient electrons are 
attracted by the charged spacecraft (ref. 5). They collide with the neutral 
atmospheric molecules In their paths and may be energetic enough to ionize the 
neutrals to form new electrons. and Ions (ref. 6). These newly created charges 
alter the space charge current arriving at the spacecraft- and shift the poten- 
tial to a lower value. The beam electrons are assumed to be energetic enough 
to leave the spacecraft completely and to play a negligible role In the Ioniza- 
tion. This mechanism Is capable of explaining the nonmonotonic current-voltage 
behavior observed. 

In this paper, we apply the sheath Ionization model to large spacecraft 
In the Ionosphere, in particular. It Is Important to find out whether the 
nonmonotonic current-voltage behavior during electron beam emissions would 
still be present for large spacecraft. Details of the method are given and 
followed by a discussion of results. 
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SYMBOLS 


e electron charge 

t electric field- 

Ijj beam current 

ra mass. of electron 

f? 

m^ mass of. Ion - 

n^ density of ambient- electrons 

n density, of Ionization Ions 

n density of Ionization electrons— — 

P probability of Ionization _ 

r radius of spacecraft 

r radial- position .measured from center of spacecraft 

r' radial position used as. Integration variable 

r 0 radius of sheath measured from center of s pacecra ft 

v e velocity of an elec ron In sheath 

v th th erroa l velocity of ambient electron 

v $ sweep velocity 

e Q permittivity of space 

X mean free path of electron neutral collision 

p space charge density 

* — electric potential 


MATHEMATICAL FORMULATION 

, rpf °LjE proach u ! ed 1s .f° st “ d y an analytical "plasma probe- model 

3 »*» 7 * 8 '? w * th space charge flow of electrons accelerating through the 
sheath surrounding a spherical "probe," which represents a spacecraft In an 
lonlzable plasma environment. Magnetic field effect Is Ignored In this model. 

niof J5 e a b 5 a ? 1S ? ss H med *° u be energetic enough to leave the spacecraft com- 
pletely and Is not stopped- by its own space charge at all. As the beam elec- 

fch«fM a J!' 4 : the spacecraft becomes charged oppositely. A polarization region 
(sheath) Is formed In the vicinity of the spacecraft. In our model-. Ions are 
assumed to be depleted due to charge repulsion Inside the sheath (fig. 1). 

The depletion radius r 0 will be defined by the balance of the outgoing 

beam current with the Incoming ambient current. For a beam current Ih the 
depletion radius r 0 Is determined by current i b , tne 


lb ■ 4«r£n e ev tb 


0 ) 


where vt b Is the thermal velocity and 
ambient electrons. Some typical values 
of equation (1) are shown In figure 2. 


n e Is the number density of 
of sheath radius as calculated by means 
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The potential 4» at any point inside the sheath is governed by 
Poisson's equation: 


where fi is the space charge density and c 0 is the permittivity of empty 
space. 


Spherical Symmetric System 

To simplify the geometry, we assume spherical symmetry in the spacecraft 
and sheath system. Equation (2) becomes simply a radial equation: 



of the potential 


4» gives the electric, field 


3<b( r 1 
ar 


E(r) 


(4) 


Talcing into account the electron and ion pairs created as a result of ioniza- 
tion, the charge density p at any point r In the sheath- is given by the 
sum of charge densities (fig. 1): 

*'r) = e (n*(r) - n~(r) - n e (r>] (51 

where n e Is the return current (primary) electron density and n f and 
n- are the ionization Ion and electron densities, respectively, due to 
return current electron collisions with neutrals. 

The ionization electron. density n~(r) Is due to all Ionizations that 
occur outward of r, and the density n + (r) of Ions at r Is due to a.11 
ionizations that occur Inward of r. Thus 


and 


n"(r) 




[2e|«>(r) - <D(r')|/m e ] 1/2 


( 6 ) 


n + (r) 



M r' 2 dr' 
ut • 

[2e|4»(r) - *(r • ) I/m, ] 1/2 


(7) 
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where 





« \lie[v e (r*)]n e (ri)v e (r*) 


( 8 ) 


Numerical -Method 

TO solve the system of equations (3) to (8), one divides the space of the 
sheath Into M . concentric shells and sets up N equations for the- N unknowns 
♦l (fl9* 4). In view of the complexity of the Ionization terms in equations 
(6) and (7), It Is Impossible to solve these equations exactly. Instead, one 
seeks the approximate solutions that minimize a function F, the mean square 
of f^, constructed from the radJal Poisson equation (eq. (3)) for the 1 th 
cell, where 1*1 N. 

)] 4r (9) 

J' 

where the electric field E (eq. (4)) is constructed In a finite difference 
scheme: 


^1 ( ^i * • • • *Ej|) * (r “ (r 




^ * • • • • Ejj 


♦i • ♦hi 


Ar(E, + 2E 


HI 


* W 


( 10 ) 


The numerical method used to solve equations (8.) to (10) Is the standard 
Newton- Raph son method of Iteration: 


E 


<M) F 0) 

1 fc 1 


•( J ) 


.E 


(J)! 


Of 4 


(J) 


E U>1 

N J 


/3E, 


(ID 


A Set of trial solutions Is used to start In the Newton-Raphson Iteration 
process, and a convergent set of solutions Is sought for each set of Input 
parameters such as beam current, ambient electron density, ambient electron 
temperature, mean free path, and spacecraft radius. 


RESULTS AND DISCUSSION 

Figure 5 shows the computed results of spacecraft potential as a function 
of electron beam current for various electron densities, electron temperatures, 
and mean free paths. The nonmonotonic behavior of potential current curves 
shows up. At low currents, the potential Increases with beam current. When 
the current Increases further, ionization occurs Inside the sheath. The 
potential then turns around -as the current of the electron beam Increases. 

The Ion and electron charges created by Ionization alter the behavior of 
the space charge flow, originally governed by the single charged Poisson equa- 
tion. The potential turns to a lower value and stays approximately constant 
as current further Increases. 
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At this stage, the potential profile as a function of radial distance 
shows locally f-lat gradient. This Vs due to Ions created Inside the sheath 
not being able to move out quickly because of their heavy masses. If a local 
Ion charge buildup forms a potential hump. Ion motion would be two ways, and 
the theory, would then break down. 


To overcome this dlffteulty, a sweep velocity v s Is added to the Ions. 
Equation (7) beco mes 


n*(r) * 

r 



dn 

dt 


2 

r‘ dr* 


4 r 


[2el4»Ct) - 4 »(r')|yjn 1 v 2 l 1/2 


( 12 ) 


It Is argued that the motion of a spacecraft relative to Its plasma envi- 
ronment can provide such a sweep velocity v s (eq. (12)). The value of 
v s Is of the order of spacecraft velocity and Is an arbitrary Input to the 
computation. However, at a higher current, a potential hump again shows up, 
and the computation falls to converge. The technique breaks down. It is con- 
jectured that two-way space charge flows should be accommodated when a poten- 
tial hump appears. 


For Increasing spacecraft radii, the nonmonotonic current-voltage behavior 
still persists (fig. 6 ). However, Increased spacecraft radius lowers the 
maximum spacecraft potential Induced by beam emission. Also, the amplitude of 
the difference between the maximum potential and the minimum (beyond the turn- 
around) diminishes. Figure 7 shows a plot of the envelope of maximum and 
minimum potentials for various spacecraft radii. 


For a given beam current I 5 (eq. (1)), the sheath surface area remains 
constant and Is unaffected by the Increase In spacecraft radius. The sheath 
thickness (defined as the sheath radius minus the spacecraft radius), however, 
diminishes. As a result, a lower spacecraft potential Is sufficient to attract 
ambient electrons, through the sheath, for the-compensatlon of electron beam 
current leaving the spacecraft. 


Beyond the turnaround point In a current-potential curve, the minimum 
potential Is limited by the minimum energy required to Ionize a neutral mole- 
cule In the atmosphere. Since such a minimum energy Is generally of the order 
of 20 eV (ref. 6 ), the minimum potential In a current-potential curve Is 
expected to approach about 20 eV asymptotically, depending on the model of 
Ionization used. For the same reason, if the maximum potential Induced by 
beam emissions is below about 20 eV, no nonmonotonic behavior Is expected. 


Figure 7 shows- the calculated envelopes of the maximum and minimum (beyond 
turnaround) potentials for various spacecraft radii In- a given ambient environ- 
ment. The amount of Ionization becomes very small as the sheath potential 
approaches the minimum Ionization potential. The amplitude of the potential 
drop beyond the turnaround also approaches the value of minimum ionization 
energy, 

There Is another critical beam current, which manifests Itself for large 
spacecraft, but not for small ones. This current Is determined by equating 
the sheath radius to the spacecraft radius. If the sheath radius Is too small, 
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s Ions ^without being charged up° U9 fhe"potei)t1a1 i: if°tho t “ col " pen J atc tlca “ Mrts- 

ssrjj isiSr 

r ?9>«»" Jf C p!^c*f t JKhSJio?*ci{! e |,S a itfSSt1«2<l* , "FiSS?o t » 1 "£ r * H, * c ‘ 1 thr<!e 

clearly. The potentlal-versus-soacorrif* !J!i< ’ p 9ure ® shows these regimes 
the smalWadlSs regime, n I s tE fSi'I? d !* mT '* re ^t1vely flat in 
occurs (l.e., this Is the realme bolLa ?k! ) n , wh1ch saturated Ionization 

voltage curve). The second regime 1° chafacterl'zed bH^h* 1 * 1 1n 3 current ~ 
potential maximum, which Is the main f£!?»™ iP zed by ^ he Presence of the 
third regime occurs when the sPa?ec?aff u nonm ?™ ton1c behavior. The 
the sheath radius (measured from the spacecraft 3 ™!!*^?**' 1ts rad1us exceeds 

The beam loses its control of the SDace™af* 3 !L C0 !!j e !^ for a 9lvert cur i”ent. 
dominates. spacecraft potential and natural charging 
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DEPLETION RADIUS (meters) 



Figure 1. - Sheath formation during beam emission. 


ELECTRON TEMPERATURE = O.l eV 
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Figure 2. -- Parametric dependence of sheath size. 
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POTENTIAL PROFILE 
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Figure 7. - Envelope of potential extrema in fig. 6. 
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Figure 8. - Nonmonotonic behavior of spacecraft potential as a function of 
spacecraft radius, for a given electron beam current. 











